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Distinctive features of Coulomb-related emissions in peripheral heavy ion collisions at 

Fermi energies 
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Light charged particles emitted at about 90° in the frame of the projectile-like fragment in semi- 
peripheral collisions of 93 Nb+ 93 Nb at 38A MeV give evidence for the simultaneous occurrence of two 
different production mechanisms. This is demonstrated by differences in the kinetic energy spectra 
and in the isotopic composition of the particles. The emission with a softer kinetic energy spectrum 
and a low N/Z ratio for the hydrogen isotopes is attributed to an evaporation process. The harder 
emission, with a much higher N/Z ratio, can be attributed to a "midvelocity" process consisting of 
a non-isotropic emission, on a short time-scale, from the surface of the projectile-like fragment. 

PACS numbers: 25.70.Lm, 25.70.Pq 



In peripheral and semiperipheral heavy ion collisions 
at Fermi energies a sizable fraction of the emitted light 
charged particles (LCPs) and intermediate mass frag- 
ments (IMFs, Z >3) are produced at "midvelocity", i.e., 
they have velocities intermediate between those of the 
projectile-like fragment (PLF) and of the target-like frag- 
ment (TLF) (see, e.g., 0, S S, 0, S S Q and references 
therein) . The midvelocity emissions represent a new, dis- 
tinctive feature of the Fermi energy domain. However, in 
spite of many efforts, a comprehensive description of its 
characteristics and evolution with bombarding energy is 
still lacking, and also the nature of its production mecha- 
nism, whether statistical or dynamical, is still a debated 
matter. 

One of the main difficulties, especially on the low en- 
ergy side of the Fermi energy region, is separating the 
midvelocity emissions from the evaporation of PLF and 
TLF. For this purpose, one usually assumes a nearly 
isotropic distribution of the evaporated LCPs in the rest 
frame of the emitting PLF (and TLF) (indeed, for mod- 
erate spins, the out-of-plane anisotropy is rather weak). 
Then, the procedure (see, e.g., 0, 01) consists in at- 
tributing the forward emissions in the PLF frame to an 
isotropic evaporation from the PLF, extrapolating it to 
all angles and finally subtracting it from the total mea- 
sured distribution to obtain the midvelocity component. 
The resulting emission pattern for midvelocity products 
is usually displayed in the (vu,v±) plane, where vii (v±) 
is the parallel (perpendicular) component of their cm. 
velocity with respect to the PLF-TLF separation axis 
(or, sometimes, the beam axis). 

Close examination of this emission pattern suggests 



that the midvelocity emissions consist not only of a broad 
distribution roughly centered at velocities intermediate 
between PLF and TLF, but also of an anisotropic "sur- 
face" emission 0, Q located along the Coulomb ridge of 
the PLF (and of the TLF as well). One can imagine dif- 
ferent mechanisms which might contribute. Early emis- 
sions — like, e.g. pre-equilibrium particles from the very 
first phases of the collision, or particles from the hot zone 
of overlapping matter during contact — are expected to 
populate mainly the central midvelocity region, possibly 
with large transverse momenta. Later emissions in the 
separation phase (or immediately after separation) are 
likely to display an increasingly strong relationship with 
just one of the reaction partners, so they are expected to 
show a kind of partial "orbiting" and to be mainly dis- 
tributed, in a non-isotropic way, on the Coulomb ridge of 
one reaction partner (usually experiments concentrate on 
the PLF), preferentially facing the other one (the TLF). 
These later emitted Coulomb-related particles may just 
be neck remnants, left behind by a dynamical multiple 
neck-rupture process or they may be produced, after 
neck rupture, by the fast non-equilibrated decay of a pos- 
sibly strongly deformed PLF, in a process resembling a 
fast oriented fission [9|, [lfj, lll| for extreme mass asym- 
metries 0, 0]- Opposite interpretations ascribe these 
later emissions to a pure statistical evaporation from the 
PLF, however perturbed by the proximity of the TLF 
(0,111). 

In this Letter we put into evidence some distinctive 
features of this "surface" midvelocity component, which 
has been measured in a common angular range simulta- 
neously with the PLF evaporation, in particular for what 
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concerns the kinetic energy spectra of LCPs and the av- 
erage isotopic composition of the Z = 1 particles. They 
may help constraining models on mid velocity processes. 

The experimental data refer to the symmetric colli- 
sion 93 Nb+ 93 Nb at 38,4 MeV^studied with the FIASCO 
setup (for details see [!, at the Superconduct- 

ing Cyclotron of the Laboratori Nazionali del Sud of 
INFN in Catania. Attention is focused on two-body semi- 
peripheral reactions, where the PLF and TLF emission 
patterns can be clearly separated. These reactions are 
selected by requiring only two heavy (Z > 10) fragments 
in the exit channel, namely the PLF and TLF (the lat- 
ter one is efficiently detected even in the most peripheral 
collisions due to the low thresholds of the gas detectors). 

For sorting purposes, we use the "Total Kinetic En- 
ergy Loss" , a parameter obtained from the kinematic co- 
incidence analysis and defined as TKEL = Ej C n m - — 
\ A^rei) where E^™ 1 - is the center-of-mass energy in the 
entrance channel, v IC \ the PLF-TLF relative velocity and 
fi the reduced mass assuming an exactly binary reaction. 
As noted in Q , at Fermi energies TKEL does not repre- 
sent any more a good estimate of the true total kinetic 
energy loss of the collision, but it is used just as an order- 
ing parameter for sorting the events in bins of increasing 
centrality. More details on the use of TKEL as an impact 
parameter estimator can be found in Appendix A of 

As the solid angle coverage of FIASCO, although large, 
is significantly smaller than 4ir, the yields of LCPs and 
IMFs are first corrected [B| for the limited geometrical 
coverage and for the low-energy identification thresholds. 
Then, using relativistic kinematics [23| . we obtain the 
emission patterns in the PLF frame (the polar axis is 
the asymptotic PLF-TLF separation axis pointing in the 
direction of motion of the PLF). One can now examine 
the distribution of the parallel and perpendicular com- 
ponents of the reduced particle momenta with respect 
to the polar axis (£A and respectively). As shown 
in Fig. 4 of the most forward part of the emission 
can be ascribed mainly to an evaporative process. On 
the contrary, at backward angles, the yields of all parti- 
cle species present a large excess, with respect to a pure 
evaporative process, which is globally ascribed to midve- 
locity emissions, with a tail of the "surface" component 
actually extending at polar angles well below 6 = 90°. 

The kinetic energy distributions of the particles in the 
PLF frame show clear distinctive features, which are best 
observed around 90° where the contribution of the early 
central emissions is strongly reduced, an appreciable con- 
tribution of "surface" midvelocity particles is still present 
and the emissions from the TLF are negligible. For ex- 
ample, the upper panels of Fig. Q] show the spectra of 
protons at TKEL = 450-550 MeV (left) and a particles 
at 550-650 MeV (right). For comparison, the lower pan- 
els show the corresponding spectra for forward emitted 
particles. All spectra have been corrected for the angle- 
dependent effects caused by the recoil of the emitter (see 
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FIG. 1: (color online) Kinetic energy spectra, in the refer- 
ence frame of the PLF, for protons (left panels) at TKEL = 
450-550 MeV and a particles (right panels) at 550-650 MeV. 
Upper and lower panels refer to particles emitted in the an- 
gular range 85° < 9 < 95° and 5° < 9 < 25°, respectively. 
The full (red) curves are fits to the data with one effusive 
Maxwell function (lower panels) and with the superposition 
of two effusive Maxwell functions (upper panels), indicated 
by the dashed (black) and dotted (blue) curves. 

[|| for more details). Such recoil effects are however rel- 
atively small at forward angles and become negligible 
around 6 = 90°. 

The distributions of forward emitted particles of the 
lower panels in Fig. [1] have been fitted with a single 
effusive Maxwell function (red curve) 

fo = N [(E - ZB )/T*] cxp[-(£ - ZB )/T ], (1) 

with Z charge of the emitted particle, Bq Coulomb bar- 
rier divided by Z and To (inverse) slope parameter of 
the distribution. On the contrary, the distributions in 
the upper panels cannot be fitted with a single Maxwell 
function, but they need a two-component function 

fu = Ni {{E - ZB x )lTl\ exp[-(E - ZB 1 )/T 1 ] 
+ N 2 [{E-ZB 2 )/T*} exp[-{E-ZB 2 )/T 2 ] (2) 

with two different slop e parameters T\ and T 2 (Ti < T 2 ) 
and equal barriers [241 ] . These two components are shown 
in Fig. [T]by the dashed (black) and dotted (blue) curves, 
respectively (the red full curve is their sum) . No attempts 
were made to reproduce the smearing of the barriers, so 
the low-energy part of the histograms (below ss 70% of 
the peak height) is not used in the fit. Fits of similar 
quality are obtained for all LCPs at all TKEL values. 

By simple inspection one sees that the slopes of the 
two components are quite different. For a quantitative 
comparison Fig. [2] presents Ti, T 2 and To for all LCPs as 
a function of TKEL. The error bars are only statistical; 



3 



14 
12 

> 10 
| 8 

H 6 

4 
2 

14 
12 

> 10 
| 8 

H 6 

4 
2 
0, 







■ □ T2 


• t o : 


p 


- 




g U 


. • ; t * 


1 : : 







I I • 



: * * 



a 



200 400 600 800 200 400 600 800 
TKEL (MeV) TKEL (MeV) 

FIG. 2: Values of the slope parameters To (filled circles) T\ 
(filled squares) and T2 (open squares) as a function of TKEL 
for protons, deuterons, tritons and a particles. The bars indi- 
cate just the statistical errors. Dashed lines through the filled 
and open squares are just to guide the eye. 



further uncertainties (due, e.g., to efficiency and recoil 
corrections and to the choice of the fit region) are around 
0.1-0.2 MeV for T and 0.5-1.0 MeV for T x and T 2 . 

The values of the slope parameters of the harder com- 
ponent, T2, are indeed definitely larger than T1 and To. 
They are similar for all Z = 1 isotopes (around 8-11 
MeV), while a still higher value (about 13 MeV) is found 
for a particles and also for IMFs (not shown in the fig- 
ure; in this latter case, a fit with a single T2 component 
is more appropriate, as the evaporation of IMFs from 
the PLF is negligible). These values compare well with 
those of [l6j], which however refer to the central part of 
the midvelocity emissions in a lighter system at higher 
bombarding energies. 

The second item worth noting is the quite similar val- 
ues of To and T\. This leads to the conclusion that the 
softer component in the spectra at 90° (slope T\ , black 
dashed curves in the upper panels of Fig. [TJ) can be at- 
tributed to the usual evaporation from the excited PLF. 
Actually, in the tail of the spectra at forward angles 
(lower panels of Fig. QJ) there is a small excess of en- 
ergetic particles (for protons and a particles, but not for 
deuterons and tritons), which might be due to a small 
pre-equilibrium contribution. A two-component fit with 
Eq. © would give a slope parameter of the evaporative 
component slightly smaller than the T obtained with 
Eq. (UJ) and hence in even better agreement with Ti, at 
the expense of an increased uncertainty in the parameters 
of the second component. 

The slope parameters of the evaporative process, To 
and Ti, display a weak, almost linear, increasing trend 
with increasing TKEL. For all particles, the values are 



rather similar and vary in the range from 3 to 6 MeV. 
These values and their dependence on the violence of 
the collision arc in good agreement with those of other 
studies focused just on the decay of projectile residues 
in peripheral collisions at Fermi energies [17l . [H| . In the 
first bins of TKEL, the values of T\ are somewhat larger 
than those expected on the base of a simple Fermi-gas 
formula, using the estimate of E* from However, T\ 
and T2 are just slope parameters. A connection with nu- 
clear "temperatures" is particularly difficult in peripheral 
collisions at Fermi energies, due to the finite width of the 
distributions both in initial mass and excitation energy of 
the emitters and to uncertainties in the determination of 
the origin of the PLF frame. These effects are discussed 
in detail in Ref. [l|| , with regard to the decay of PLFs 
produced in peripheral collisions at Fermi energies and 
detected with the INDRA multidetector. 

What we want to stress is the remarkable difference be- 
tween the slopes T\ and T2, which is found for all LCPs 
irrespective of the value (and hence of the selectivity) of 
the adopted ordering parameter. This fact indicates that 
the physical mechanisms associated with the two com- 
ponents always maintain different characteristics. Going 
toward smaller values of the angle 8 in the PLF frame, 
T\ does not change, while T2 shows a weak tendency to 
decrease, but it remains definitely higher than T\ . At the 
same time, the percentage of Coulomb-related midveloc- 
ity emission decreases. Around 90° it amounts to about 
25, 50, 60 and 30% of the total emissions of p, d, t and 
a, respectively, and it decreases to less than 10% in all 
cases for 6 < 40°. 

If the two components observed in the kinetic energy 
spectra correspond indeed to different production mecha- 
nisms, one may expect them to differ also in other proper- 
ties, like in the average isotopic composition of the emit- 
ted particles. For this purpose we have estimated the 
N/Z ratio of the hydrogen isotopes (for which the detec- 
tors have good isotopic resolution) in the angular range 
85° < 9 < 95° in the PLF frame. One can proceed in two 
different ways. The first one is similar to that adopted 
in Ref. One first builds the whole evaporative in- 

variant cross section, cr^ p (-"-, ^), by rotating the part 
measured in the forward direction so as to fill the whole 
angular range around the PLF position. Then, one sub- 
tracts the so determined of™ p from the measured total 
invariant cross section of particles emitted in the forward 
cm. hemisphere, af°', thus finally obtaining an estimate 
of the yields of midvelocity emissions of all types. The 
results are shown by the open triangles in Fig.[3j The sec- 
ond completely different way to proceed is to deduce the 
average N/Z of the midvelocity emissions directly from 
the fitted intensity parameter N%: the results are shown 
by the open squares in Fig. [3l This second method gives 
results in agreement with the previous one, although with 
larger uncertainties. The other fit parameter, N±, is used 
for an estimate of the isotopic composition of the evap- 
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FIG. 3: N/Z ratio for hydrogens as a function of TKEL. Open 
and full squares refer to the midvelocity and evaporative com- 
ponents, respectively, and are obtained from fits of the kinetic 
energy spectra at 85° <6< 95° in the PLF frame. Open tri- 
angles refer to the midvelocity component in the same angular 
range and are obtained by subtracting the evaporative com- 
ponent (see text); the dashed line is a quadratic fit to guide 
the eye. Full and open circles (with full lines) show the results 
of 0] for evaporation and total midvelocity, respectively. 



emissions (with a clear origin from the PLF "source" , as 
manifested by their concentration on the Coulomb ridge 
[H, [f|) consist of two components. The softer compo- 
nent can be identified with the usual evaporation from 
an equilibrated source, on the basis of both the particle 
kinetic energies and the N/Z ratio of the hydrogen iso- 
topes. The other component, which displays harder ki- 
netic energy spectra of the particles, shows a more exotic 
N/Z, with a clear neutron enrichment for bound neutrons 
with respect to evaporation. These results represent a 
benchmark against which models describing midvelocity 
processes should be tested. In any case the presented re- 
sults are compatible with a picture in which the midve- 
locity "surface" emission comes from the highly excited 
contact region between PLF and TLF. After the neck 
rupture, the PLF may remain partially deformed and lo- 
cally highly excited; as a consequence, it may tend to 
evaporate particles and fragments with high kinetic en- 
ergy and with an isotopic composition resembling that of 
the neck region, before reaching full equilibration. 



orative component, shown by the filled squares. Finally, 
in the same figure, the filled and open circles show the 
N/Z values obtained in Ref. [5| for the evaporative and 
total midvelocity emissions, respectively. 

The average N/Z for the soft component of the ki- 
netic energy spectra (filled squares) nicely agrees with 
the value obtained for the evaporation at forward angles 
(filled circles), thus landing further support to the iden- 
tification of the soft component (with slope parameter 
Ti) with an evaporative process. As shown in [f|, the 
experimental values are indeed in close agreement with 
calculations employing the statistical code Gemini [20| . 

For the Coulomb-related part of the midvelocity emis- 
sions, the open triangles are in a position intermediate 
between the open and filled circles at all TKEL. This in- 
dicates that this Coulomb-related part of the midveloc- 
ity emissions, although not as neutron rich as the bulk of 
the midvelocity emissions from the central "source" , nev- 
ertheless presents an average isotopic composition defi- 
nitely higher than that of the usual evaporation (at least 
for the Z=l particles). This fact, while confirming a sub- 
stantial difference between the two mechanisms, may be 
taken as an indication of a possible enrichment in bound 
neutrons of "midvelocity" matter [E S 0, 0, EH [12], 
and/or it might be related to the reduced size of this 
"source" and hence to its higher energy concentration 

1- 

In conclusion, we have investigated some aspects of 
LCPs emitted in peripheral collisions at Fermi energies. 
We have shown that, still at 90° in the PLF frame, these 
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